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Abstract 

Using a suitable change of probability measure, we obtain a novel Poisson series representation for the arbitrage- 
free price process of vulnerable contingent claims in a regime-switching market driven by an underlying continuous- 
time Markov process. As a result of this representation, along with a short-time asymptotic expansion of the claim's 
price process, we develop an efficient method for pricing claims whose payoffs may depend on the full path of the 
underlying Markov chain. The proposed approach is applied to price not only simple European claims such as 
defaultable bonds, but also a new type of path-dependent claims that we term self-decomposable, as well as the 
important class of vulnerable call and put options on a stock. We provide a detailed error analysis and illustrate the 
accuracy and computational complexity of our method on several market traded instruments, such as defaultable 
bond prices, barrier options, and vulnerable call options. Using again our Poisson series representation, we show 
differentiability in time of the pre-default price function of European vulnerable claims, which enables us to rigorously 
deduce Feynman-Kac representations for the pre-default pricing function and new semimartingale representations 
for the price process of the vulnerable claim under both risk-neutral and objective probability measures. 
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1 Introduction 



Regime-switching models are aimed at capturing the appeaUng idea that the macro-economy is subjected to regular, yet 
unpredictable in time, regimes, which in turn affect the prices of financial securities. For example, structural changes 
of macro-economic conditions, such as inflation and recession, may induce changes in the stock returns or in the term 
structure of interest rates and, similarly, periods of high market turbulence and liquidity crunches may increase the 
default risk of financial institutions. This has been empirically verified in the stock market, as stated by |Ang and 



Bekaert ( 2002-b I , who proposed the existence of two regimes characterized by different levels of volatility. Similar 



findings have also been documented in the bond market (see Ang and Bekaert (2002-al, Ang and Bekaert (2002-c 



and Dai et al. (2007)). Most recently, in the credit market, Giesecke et al. (20111 suggested the existence of three 
regimes, associated with high, middle, and low default risk, via an empirical analysis of the corporate bond market 
over the course of the last 150 years. These considerations have led many researchers to use regime-switching models 



for asset pricing. In the context of options on stocks, Guo (20011 considered a market consisting of two regimes, and 



provided a semi-analytical formula for the option price, based on occupation time densities. |Buffington and Elliott| 
(2002) generalized the method to the case of multiple regimes, under the assumption that the generator of the Markov 



chain is time homogenous, and derived expressions for the price of European claims. Gapeev and Jeanblanc (2010) 
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obtained closed form expressions for European claims, assuming geometric Brownian motion dynamics with regime- 
switching drifts, while Yao et al. (2006) developed a fixed point iteration scheme to recover prices of European options 
in a regime-switching model, assuming time homogenous generators. Regime-switching models for default-free interest 
rates derivatives have been studied by EUiott and Wilson^ (2007| , Elliott and Siu (2008), and Siu (2010-a). |Norberg| 
( 2000 ) studied no-arbitrage pricing of derivatives on the regime parameters with emphasis on computation. 

Most of the literature has focused on default-free European style claims with "simple" payoffs (i.e. payoffs depending 
only on the value of a random variable such as the terminal value of the underlying asset or the underlying source of 
randomness) with the exception of a few works detailed next. Guo and Zhang ([2004^) studied American- type default-free 
claims under regime-switching models. Graziano and Rogers (2006) provided a methodology to price barrier options 
with a regime-switching dividend process. In the context of credit risk, Siu (2010-b) considered pricing of interest rate 
derivatives with a defaultable counterparty. Siu et al. (2008) analyzed credit default swaps, and derived a system of 
coupled partial differential equations yielding the term structure of default probabilities needed for pricing. [Bielecki] 



et al. (2008-al and Bielecki et al. (2008-b) studied pricing and hedging of a defaultable game option under a Markov 



modulated default intensity framework. Blanchet-Scalliet and Jeanblanc ( 2004 ) gave the hedging strategy of vulnerable 
claims using defaultable zero-coupon bonds and default-free assets, under the assumption that the default-free market 
is complete. Bielecki and Jeanblanc (2005) applied the theory of indifference pricing for pricing and replication of 
vulnerable claims within the reduced form framework. For an excellent survey on pricing and hedging aspects of 
vulnerable claims, the reader is referred to Bielecki et al. (2004). 

In this paper, we consider a macro-economy with finitely many observable economic regimes, containing state 
information regarding the equity (drift and volatility), credit (hazard intensity and loss given default), and interest 
rate market (short rate). We consider three liquidly traded securities, namely, a money market account, a risky 
(default-free) stock, and a defaultable bond. The dynamics of the securities are assumed to depend on the macro- 
economic regimes, which are modeled using a finite state continuous time Markov process. We follow the reduced form 
approach to credit risk and model the default event using a doubly stochastic framework. Pricing of contingent claims 
is still problematic in the context of regime-switching models, in part, due to the lack of easily computable expressions 
for option prices. Hence, explicit characterizations of these quantities that allow their efficient evaluation are of great 
need in the field. Our paper contributes to this fundamental problem by developing a novel approach for pricing 
vulnerable claims in a regime- switching model consisting of an arbitrary finite number of regimes, whose dynamics is 
governed by a possibly time varying generator. The proposed method can also be applied to price a certain type of path 
dependent claims, referred hereafter as self-decomposable claims, whose payoffs can be decomposed in terms of payoffs 
of shorter maturity claims, of possibly different type (see Section 5.1 for the precise statement). Self-decomposable 
claims encompass not only basic instruments, such as bonds, whose price may also be recovered from no-arbitrage 
arguments via the solution of a coupled system of ordinary differential equations (ODE's), but also more exotic path- 
dependent instruments, where prices can only be recovered via Monte-Carlo methods. As an example, we apply our 
algorithm to price a type of barrier options on the volatility of the stock, which turn out to be self-decomposable in 
terms of shorter maturity barrier options and bond prices. Note that, in our macro-economic model, our regime driven 
parameters for the short rate, volatility, and the default intensity may be seen as "quantized" proxies for the actual 
interest rates, volatilities, and credit spread indices, and, therefore, options on these underlying financial measurements 
provide means for the investor to hedge against interest rate risk, market and default risk, in different economic regimes. 

The key tool behind our results is a type of Poisson series representation for the prices of vulnerable claims obtained 
by applying a novel change of probability measure technique, which transforms the underlying Markov process into a 
"homogenous" process with constant transition intensities. This Poisson representation enables us to express the price 
of the claim as a series expansion in terms of Laplace transforms of the symmetric Dirichlet distribution, which can be 
efficiently evaluated by Taylor approximations whenever the claim's maturity is small. Hence, a second important step 
in our method is to decompose the claim's price in terms of shorter dated claims. We provide a unified framework to 
quantify (and hence to control) the error due to the Taylor approximations, the error due to the discrete approximation 
of the time-varying generators, and even the error due to the truncation of the Poisson series representation. Our 
algorithm is computationally fast and, even for simple claims such as bonds, is able to achieve a high level of accuracy 
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in the price, within a time complexity which compares favorably with standard ODE methods. As a far reaching 
application of our method, we also propose a new method to price vulnerable call/put options on the stock (where 
there is an additional risk factor represented by the Brownian motion), a subject which has received significant attention 
in the literature, as documented above. 

Under mild assumptions on the Markov chain generator, we also apply our Poisson series representation to give a 
rigorous proof of the differentiability of the pre-default price function. This, in turn, allows us to provide Feynman-Kac 
and semimartingale representations for the pre-default price process (given by Eq. ([9]) below), and for the vulnerable 
claim price process (given by Eq. (10) below). This extends earlier results in the literature which, on one hand, had 
focused mostly on martingale representations and, on the other, typically presumed the time differentiability of the price 
function (e.g., Elliott and Kopp (2004), Elliott et al. (2007)). Needless to say that our semimartingale representations 
play a key role in the hedging of defaultable contingent claims, a problem which has attracted increasing attention 
since the start of the credit crisis, due to the high number of credit quality deteriorations, or even default events, 
experienced by financial institutions. Although we do not explicitly address the hedging problem in this paper, our 
formulas can be used within any hedging framework employing Markov-modulated dynamics. We refer the reader 
to Bielecki et al. 



and 



(2003-a), Bielecki et al. (2003-b), Bielecki et al. (2004), Blanchet-Scalhet and Jeanblanc (2004), 



Bielecki and Crepey ( |2011 ) for hedging methodologies based on the mean- variance or BSDE approach. The 



semimartingale representations are also fundamental in portfolio optimization problems as demonstrated in [Capponi] 



and Figueroa-Lopez (2012). Another related work is Liew and Siu (2011), where a martingale representation for a 



two-dimensional martingale generated by the Brownian motion and the Markov chain was provided with the objective 
of establishing the existence of an admissible investment strategy in a default-free regime-switching market. 

The rest of the paper is organized as follows. Section [2] sets up the defaultable regime-switching model. Section 
|3] establishes the Poisson series representation of the vulnerable claim price function. Section |4] provides rigorous 
semimartingale representations for the price process of the vulnerable claims. Section [5] presents our novel efficient 
method for pricing vulnerable claims and provides a rigorous error analysis. Section [6] illustrates the accuracy and 
computational power of our algorithms on concrete choices of claims. Section [7] concludes the paper. Proofs and 
numerical details are deferred to the Appendix. 



2 The defaultable regime-switching model 



As in Capponi and Figueroa-Lopez (2012), we consider a regime-switching market model with a defaultable security 
driven by a continuous-time Markov process X. More specifically, our market consists of a risk-free asset, a risky 
(default-free) asset, and a defaultable T-claim (such as a bond) with respective price processes {-£^£}i>05 {*^t}i>0: and 
{7(t, r)}o<t<T defined on a complete filtered probability space (17, Q, G, P). Here, P denotes the real world or historical 
probability measure and G := [Qt) is an enlarged filtration given by Qt '■— J-t V Ht, where F := {Ft\t>o models the 
reference filtration and T-Lt ~ <j{H{u) : u <t) \s the filtration generated by an exogenous default process H{t) := lr<t, 



after completion and regularization on the right (see Belanger et al. (2004) for details). Throughout, t represents the 



default time of the defaultable security of interest. We also adopt the canonical construction of the default time r in 
terms of a given hazard process {/it}t>o, so that 



inf U e R I 



h,,du > Q 



(1) 



where is an exponential random variable defined on the probability space (ri,t/,P) and independent of F. In that 
case, it follows that 



(.J:^H{t)- [ {l-H{u-))Kdu 
Jo 



(2) 



is a G-martingale under P (see Bielecki and Rutkowski (2001), Section 6.5). The intuition of ^ is that one needs to 
compensate the single jump process for default prior to the occurrence of a default. Let us also recall that, under the 
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previous setting, F satisfies the so-called martingale invariance property with respect to G, which states that every 



F martingale is also a G martingale (see Sections 8.3.1 and 8.6.1 in Bielecki and Rutkowski (2001)). The previous 
principle is also typically referred to as the H-hypothesis . 

As previously mentioned, we place ourselves in a regime-switching market model. More specifically, we define 
an F-adapted continuous-time Markov process {Xt}t>a with finite state space {ei,e2, . . . ,eAr}, where hereafter ei = 
(0, 1, ...0) e and ' denotes the transpose. Throughout, Pij{t, s) :— V{Xs — ej\Xt — Ci) ior t < s represents the 
transition probabilities of X and A{t) [ai,j(i)]ij=i,2,...,Ar denotes the generator, defined by 



it) = lim 



Pi^j{t,t + h) 



ai,i{t) := -^aij(i). 



The following scmimartingale representation is well-known (cf. Elliott et al. (1994)) 



Xt^Xa+ I A'{s)Xsds + Af{t), 
Jo 



(3) 



(4) 



where M (t) — (Mf (t), . . . , M]i^(i))' is a R -valued P-martingale process. The following terminology will also be 
frequently needed in the sequel: 



N 



(5) 



We now specify the three instruments in the market, whose dynamics are driven by {Xt}t>o- We have a locally 
risk-free asset {Bt}t>o with Markov-modulated dynamics 



dBt = rtBtdt, 

where rt takes a constant value r-i if the economy regime variable Xt is at the i*^ state e^. That is. 



(6) 



rt 



{r-Xt) 



Ct ' 

where throughout (•,•) denotes the standard inner product in M.^ and r = (ri, r2, . . . , r^y)' is a vector of positive 
constants, which represent the possible risk-free short-rates in the market. The risky (default-free) asset {>5't}t>o 
follows the dynamics 

dSt ^ fj-tStdt + (JtStdWt, Sq^s, (7) 
where {VFt}t>o is an F-adapted standard Brownian motion independent of {Xt} under the measure P, and 



Ht (/i, Xt) = ■- 



c^t 



(8) 



for constant vectors /i — (/ii, /i2, . . . , Mw)' ^^'^ — ("'ii '^2, • ■ • , o'nY , representing the respective appreciation rates and 
volatilities that the risky asset can take depending on the different economic regimes. 

The pricing methods we introduce in this paper can be applied to a wide range of defaultable securities such as 
corporate bonds, vulnerable claims, recovery payments at default, or even vulnerable call/put options on the risky 
stock {St}t>o- As a way to fix ideas, we will focus on a defaultable security, whose recovery payment at default is 



determined by the recovery of market value (RM) assumption (see, e.g., (McNeil et al. 2006 Section 9.4) or (Bielecki 



and Rutkowski 2001, Section 8.3)). More precisely, this security delivers an J"T-iiieasurable promised payment X at 
time T if no default has occured by that time (i.e. t > T) and, otherwise, delivers the recovery payment (1 — Lt-)7(t^; T) 
at the default time t < T, where {jit; T)}t>o is the pre-default value process of the security and {it}t>o is an F-adapted 
process, which represents the percentage loss given default. Then, the following fundamental risk-neutral valuation 



formula due t0 |Duffie and Singleton ( 1999 ) (see also Proposition 8.3.3 in [Bielecki and Rutkowski ( 2001 )) is well-known: 



;Y(,-St(r,+h,L,)ds 



J't 



(9) 
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where Q is the equivalent martingale measure used in pricing. Furthermore, denoting {T{t;T)}t>o the value of the 
previously described defaultable security, it follows that 



An important example of the defaultable security (|10|) is a defaultable bond given by the price process 

g-/,^(r.+/i,L,)ds 



p(t,T) = l{,>4}E^ 



(10) 



(11) 



Throughout, we assume that Q is such that W is still a standard Brownian motion and X is an independent 
continuous-time Markov process, under Q, with possibly different generator :— {of ^(t)\i.j=\ ■ The existence 

of the measure Q follows from the theory of change of measures for continuous-time Markov chains (see Section 11.2 
in Bielecki and Rutkowski (2001 )). Such density transformation theorems for Markov chains will play a key role in our 
pricing approach below. Concretely, for any i ^ j ^-nd some bounded measurable functions Ki.j : M+ — > (—1, oo), let 



:=a,j(t)(l + «,j(i)), 



(12) 



and, for i ~ j, let 



AT 



We also fix Ki^i(t) = for i = 1, . . . , A^. Now, consider the processes 

ft 



Hi' 







aij{u)Hldu, 



where 



:= l{Xt=e.}, and -ff^'^ ^ l{jf„_=e.}l{x„=e,}, (i^j)- 



(13) 
(14) 



0<u<t 



The process {Ml'-'}t>o is known to be an F-martingale for any i j (see Lemma 11.2.3 in 



Bielecki and Rutkowski 



(2001)) and, since the _ff-hypothesis holds in our default framework, they are also G-martingales. Then, we have the 



following result (see Proposition 11.2.3 in Bielecki and Rutkowski (2001)) 



Lemma 2.1. The probability measure Q on G = {Gt)t with Radon- Nikodym density {r]t}t>o given by 



1+ / E ^--^-j(")^^^«'' ^^0' 

Jio.t] ,^^.^1 



(15) 



is such that X is a Markov process under Q with generator [a^.(t)]ij=i,2,...,Ar- 



Without loss of generality, one can take the measure Q in the previous lemma to be such that W is still a Wiener 
process independent of X under Q. Analogously to (|4|, the process 



M^{t) -.^Xt- Xo 



syXsds, 







is a 11 -valued martingale under Q. Comparing (16 1 to (|4|), note that 

~-M^{t)+ f {A{sy - A^{sy)Xsds. 
Jo 



(16) 



(17) 
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To be consistent with our regime-switching market model, we assume that the default intensity h and loss rates L 
are also driven by the underlying Markov process X] i.e., 

ht := {h, Xt) - h^^ , and Lt (L, Xt) = L^^ , (18) 

for constant vectors h = {hi, ... , hfq)' and L = {Li, . . . , i^v)', where hi > 0, and < < 1. We emphasize that the 
distribution of the hazard rate process ht = {h,Xt) under the risk-neutral measure is different from that under the 
historical measure, which allows to incorporate a default risk premium into our framework. 

The following result states that the process {^f }t>o introduced in ^ is also a Q-martingale. However, in order to 
indicate in the sequel when certain dynamics are being taken under Q or under P, we will introduce a new notation 
{C?}t>o- We should keep in mind throughout that = The proof of this result is presented in Appendix [a] 

Lemma 2.2. The process 

if -.^Hit)^ [ {l-H{u-))h^du (19) 
Jo 

is also a (local) martingale. 



3 Poisson series representations for vulnerable claims 

In this section, we introduce the main tool upon which the construction of our novel algorithms as well as the semi- 
martingales and Feynman-Kac representations of the price processes are based on. We focus on the case when the 
promised payoff X in ( 10 ) is of the form S(Xt) for a deterministic function S : {ei, 62, ... , ejv} — > K. In this case, due 



to the Markov property, the pre-default value process (|9| can be written as 
where, for each economic regime i G {1, . . . , N}, 



(20) 



The key idea of our approach to compute ( 20 1 lies in changing the risk-neutral probability measure Q into an 
equivalent measure Q such that {Xt}t<T becomes a homogeneous Markov process under Q. Such a probability measure 
Q exists whenever afj{t) is strictly positive for alH > and i 7^ j. Concretely, define functions kij : [0, 00) (—1, 00) 
such that 



1 



7V- 



for any i ^ j, and 



= 0. 



IS a 



We also let Uij :— 1/{N — 1) for any i ^ j and = —1 for any i — 1,...,N, so that A :— [d 
valid generator of a homogeneous Markov process with transition times determined by a homogeneous Poisson process 
{Mt}t>o with intensity 1 and an embedded Markov chain {Xi}i>i with transition probabilities pij := 1/{N — 1) for 
i ^ j. Now, let us define a probability measure Q with Radon-Nikodym density {77t}t>o given by 



Vt 



N 

Vu-K^,J{u)dM'^^, 



io.t] 



(21) 



where M^'^ := Hl'^ 



{u)Hl^du, and we used notation (14 1. By virtue of Proposition 11.2.3 in 



Rutkowski 
follows: 



/o ^ij > 

(2001), {Xt}t>o is a continuous Markov process with generator A under ' 



Bielecki and 



The representation below then 
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Theorem 3.1. Suppose that, for any i 7^ j , the function af^ defined in (12) is such that 



0< inf |a^^.(.)|< sup |a^,-(.)|<oo. 



The 



Xt = e,: 



(22) 



(23) 



(24) 



As explained before, Q has the virtue that the generator of {Xt]t>G under Q is given by 5;^- — \/{N — 1) {i ^ j) 
and i = —1. Thus, by conditioning on the number of transitions in {t,T] (which is necessarily Poisson distributed 
with unit intensity), we get the following series representation: 



where K{t) :— [Ki,j{t)]ij and r{t) :— (fi(t), . . . , fjv(t))' are defined as 

K^Jt) := - log ({N - l)af^{t)) n{t) := n + h,L, - 1 - af/t) 



where, denoting [•] = 



\Xo = e, 



-(^-*)^^^^$.™[S](r-t;T), 



m=0 



ml 



(25) 



S(X„) expl^J2 r{T~C + s)' X^^ds - ^ X'^_^K{T - ( + CU^„))Xn 

[ n=0-^?'^(") n=l 



(26) 



Here, {Xi} is the embedded Markov chain of {Xt}t>o and C/(i) < C/(2) < • • • < C/(r,i) are the ordered statistics of m 
i.i.d. uniform [0, 1] variables independent of {Xi}, fixing C/(o) •= ^-i^d C^(m+i) •= A\ Thus, for instance, we have: 



$,,o[S](C;T)=S(e,), $^,i[2](C;T) 



N - 



(27) 



T-C 



The series representation (25) is one of the fundamental tools in the forthcoming results. We shall apply it to find 



conditions for the pre-default price functions (20) to be differentiable, which in turn will lead to rigorous Feynman- 



Kac and semimartingale representations for the pre-default price functions ^ and the vulnerable claim price process 
(10), respectively (see Section [4] below). Furthermore, we shall also apply the above Poisson series representation in 
developing our pricing methods for vulnerable contingent claims (see Section [5] below). In what follows, the dependence 
on T will sometimes be omitted from $ in order to lighten the notation. 



4 Semimartingale representations of vulnerable claims 



In this section, we rigorously deduce Feynman-Kac martingale representations for the pre-default price process ( 20 ) , and 



semimartingale representations for the price process of the vulnerable claim given by Eq. ( 10 ) . The latter are given both 
under the risk-neutral probability measure Q and under the historical probability measure P. The representations will 
be expressed in terms of the martingales ^ and (|4| (respectively, (16) and (19)). The semimartingale representations 
obtained here are vital in many fundamental problems of finance including hedging and portfolio optimization. Hedging 



with credit risky securities has attracted a lot of attention in recent years (see, e.g., Bielecki et al. (2003-a), Bielecki 



et al. (2003-b), Bielecki et al. (2004), and Blanchet-Scalliet and Jeanblanc] (2004)), while, combined with our portfolio 

^Note that, strictly speaking, the distribution of the variables = f/(o) < ^(1) < • ■ ■ < f^(m) < f^(m+i) = 1 depends on m and, hence, 
we should probably denote them as {C/^'™'}™.]^, but for simplicity, we will omit the superscript (m). 
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optimization framework developed in Capponi and Figueroa-Lopez (2012), our results here can also be applied to 



portfolio optimization problems with more general vulnerable claims of the form ( 10 1 



The representations obtained in this section are based on the following result, which in turn follows from the Poisson 
series representation ( 25 ) introduced in the previous section. Its proof is presented in Appendix [A} 



Lemma 4.1. In addition to the conditions in 
(0, T) and that 



22), suppose that the functions af^ are continuously differentiable in 



sup 

sG[0,T] 



d4M 



ds 



< oo, for all i ^ j. (28) 
Then, the pre-default price function ^'i(t) :— ^i\B\{t\ T) defined in \2(f^ is continuously differentiable for any t £ (0, T). 



4.1 Markov chain dependent claims 

We proceed to establish a system of ordinary differential equations for the pre-default price functions ^ and the 
semimartingale representation for the price process ( 10 ) in the case of a claim X depending only on the underlying 



Markov chain. An important example of a Markov chain dependent claim is the defaultable bond price, which amounts 
to taking S(ei) = 1 for all i — 1, ... ,7V. Similar Feynman-Kac representations have been given in Elliott and Kopp 



(2004), Corollary 9.8.4. However, these are given under the assumption that the price of the regime conditioned claim is 



a time differentiable function. Here, we fully justify the differentiability assumption using our Theorem |3.1| and Lemma 
|4.1| and, thus, make the analysis fully rigorous. To the best of our knowledge, the semimartingales representations are 
new and, as previously indicated, are crucial in, e.g., hedging and portfolio optimization problems. 
The following result, whose proof is presented in Appendix \K\ is the main result of this section. 



Theorem 4.2. Suppose that the conditions (22) and (T^M) are satisfied. Then, the pre-default price function '^iit) 



'^i[^]{t;T) is differentiable for any t G (0,T). Furthermore, the following statements hold: 

(1) The vector of pre-default prices '^{t) — (^'i(i), . . . , 5'7v(i))' satisfies the following coupled system of backward 
Ordinary Differential Equations (ODE): 



*,(T)=S(e,), 



.N. 



(29) 



(2) The prices process (10) of the vulnerable claim with payoff 'B{Xt) has the following semimartingale representation: 

0<t<T, t<T, (30) 



under the risk-neutral probability measure 



(3) The prices process (10) of the vulnerable claim with payoff 'E,{Xt) has the following semimartingale representation: 



dT{t,T) ^ T{t-,T) (^[rt + ht{Lt - 1) + I?(t)] dt + ^"^^^j^j^^^^^^^^ - d^!^ . 0<t<T, t < t, (31) 
under the historical probability measure ¥, where D{t) := {(Di{t), . . . , DN{t))' , Xt) with 



N 



D,{t) 5](a,,,(<) - a%(t))^^ = Y.{a,,(t) - a%{t)) 



- 1 



(32) 
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Remark 4.3. As a direct consequence of the previous result, in the time-invariant case (i.e., a^j(t) = afj), the solution 
of (29) can he expressed in "closed- form" as 

^(t;r) =e-(^-*)^^H, (33) 
where H = (S(ei), . . . , S(eAr))' and the components of the matrix are 

[Fi,]is = ri + h^Li - -, and [F^]i,j = -af,, i 7^ j- 



4.2 Markov chain and diffusion dependent claims 



In this section, we provide the Fcynman-Kac type partial differential equations and the semimartingale representa- 
tions of the price process for claims whose payoffs depend both on the Brownian paths as well as on the underlying 
Markov chain. For non- vulnerable claims, Elliott et al. (2007) and Yao et al. (2006) obtained similar Feynman-Kac 
representations. Here, the analysis is extended to deal with vulnerable claims. 



We follow Elliott et al. ( 2005 ) and choose the the risk- neutral measure Q such that 

2 ~1 



exp ■ 







1 







ds 



In particular, the (J^t,Q) Brownian motion is given by 



rs - f-s 



ds. 



(34) 



(35) 



Note that, without loss of generality, we can take Q so that under it, {Xt}t>Q is still a continuous-time Markov process 
with generator A'^, independent of {W^(^}f>o. 

Next, we consider a pre-default price process of the form ([9]), where X — p(S't), for a deterministic function g. We 
denote it by U{t] s) = (ni{t; s), . . . , IlN{t; s))' , where 



Ui{t;s) := E*^ 



The proof is presented in Appendix section. 



Theorem 4.4. Assume that each Ili(t;s) is continuously differentiable in t ^ (f^^T) and twice continuously differen- 
tiahle in s ^ (0,cx)). Then, the following statements hold: 

(1) The vector of pre-default prices n(t; s) = (ni(t; s), . . . ,njv(t; s)) ' satisfies the following coupled system of back- 
ward Partial Differential Equations (PDE), in the domain 0<t<T,0<s< oo: 



5n,;(t;s) 
dt 



dll,{t;s) 



ds 



li{T,s) = g{s), 



1, 



(36) 



(2) The prices process (10), denoted hereafter by T{t,T; St) , of the vulnerable claim with payoff X = (?(S't) has the 



following semimartingale representation, for < t < T, < St < oo, t < t: 

dT{t, T; St) = V{t~,T- St^) [[rt + ht{Lt - l)]dt 
under the risk-neutral probability measure Q. 



(n{t;St),dM'i{t)) d{nifA).x,) ^, , , 

' '^tSt^^^yr^^^r^^dW^ - d^Y ] (37) 



{nit;St),Xt-] 



'■{n{t;St),X^^) 
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(3) The prices process ( 10) of the vulnerable claim with payoff g^Sx) has the following semimartingale representation, 



forO<t<T,0<St<oo,t<T: 

dT{t,T;St) = Tit-, T: St-) 

(Hit; St), dM^ it)) 



g 0(n(t;St),Xt) 



dt 



+ (JtSf 



d{n{t-St),Xt) 

ds 



-dWt-d^f , 



m;St),Xt-) '{n{t;St),Xt-] 
under the historical probability measure ¥, where D{t, Sf) '■= {{Di{t, St), ■ ■ ■ ,Df^{t, St))' ,Xt) with 



(38) 



(39) 



5 Novel algorithm for pricing vulnerable claims 

This section introduces a novel algorithm for pricing vulnerable claims. Concretely, Section |5.1| presents the basic 



algorithm and analyzes its computational complexity. Section 5.2 develops an error analysis of the algorithm by 
providing precise bounds for approximation error of the method. Section |5.3| then builds on the basic algorithm and 
extends it to the pricing self-decomposable claims on the underlying chain. Section |5.4| introduces a novel algorithm 
to price vulnerable European type options using the basic algorithm. 

5.1 Markov chain dependent claims 

In this section, we develop a novel efficient algorithm for pricing a claim whose payoff depends on the underlying 
economic regime in place {Xt}t>o- This algorithm is used to price simple European claims of the form 'E{Xt) for a 
deterministic function S as in Section [s] Let us consider the pre-default price function ( 20 1 under the i*'' regime. Our 
approach to compute (20) builds on the series representation (25). Concretely, we will derive a formula for $i_m[S](C) 



when the risk-neutral generator is time- invariant. The formula will be expressed in terms of the Laplace transform 
of the "symmetric" Dirichlet distribution, defined by 



(40) 



where Tm '■— {{M, ■ ■ ■ , ^m) G IR™ ■ > 0, X^I^i — The proof of the following result is given in Appendix [A| 

Proposition 5.1. Suppose that cifj{t) = afj (hence, K and r are also time-invariant). Then, for m > 1, we have 
that 



1 



{N - l)^ 



such that Cj € {ei,...,eAr} and 



S(e™)e-^^''^"-^-i'^"-i^^^"£™(Cf'(e~o-g™),...,Cr'(g„»-i-g„»)), (41) 

(ei,...,e,„) 

where we set cq = Ci and the above summation is over all "paths" (ei,...,e,i 
ij ^ ej-i, for i 1, . . . ,m. 

Remark 5.2. Note that 
<t,,„|H](0 = H(e,)e-<'', 



^ ^ je{i,...,7V}\{»} ^ 



1 - e-^^^^-^i' 



(42) 



For a general m, the following Taylor approximation around the origin will turn out to be quite useful to compute the 
I transf 



Dirichlet Laplace transform {40): 

m + 1 (m + l(m + 2 \ 



3/2 



(43) 
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In order to evaluate the option price, the infinite series (251 wih be truncated and the Laplace transform in (41) 



may also be approximated by ( 43 ) . Both of these approximations are valid when timc-to-maturity ( = T — t is small. 



Hence, it useful to express the option price (23) in terms of the price of options with shorter expiration. Concretely, 



fix a small mesh parameter S := T/k for a given positive integer k and let 

se(«/u]:AXs^O 



(44) 



Then, using the tower and Markov properties, the time price of the vulnerable claim may be computed as follows 



*4S](0;r) = [e'^>-^E{XT)\Xo = e,] = 



Xn 



^{Xs] 



Xo 



*4S](0;<5), 



(45) 



where we introduced a new payoff function S defined as S(ej) = 5'j[S](0;T — 6). Note that the last expression in 



(45) repr esen ts the price of an option with short maturity S and, hence, it can be accurately computed by taking M 
terms in (25) and (possibly) using (43). In order to evaluate the option's payoff S, we apply again the procedure (45), 
replacing T by T — (5. Computationally, one can create a recursive or iterative implementation in order to evaluate 
( 45 ) . The pseudo-code of the proposed iterative algorithm is given in Appendix [b] (Algorithm [I] therein) . 



Remark 5.3. The previous method can be adapted to deal with sufficiently smooth time dependent functions f{t) and 
Kit) (e.g., continuously differentiahle in [0,T]j. Indeed, as in (45) and using notation we have 



^^mt■.T) = EQ [e'^-^^{XT)\Xt = e,] = E^ [e^"+^^*c,+. [S](t + 6;T)\Xt = e,] - *.[S](i;< + <5), 



(46) 



fixing S(ei) :— '^i\E]{t + 6;T). Then, for sufficiently small 5, ^i\S\{t\T) can he accurately computed as if f and K 



were time-invariant during the period [t,t -\- 5]. We analyze the error of this approach in the section 5.2.2 below (see 



also the end of Section 5.2.1 for another approach) 



Remark 5.4. The complexity of the proposed method (Algorithm^in Appendix\B^ can be evaluated as follows. Let 



$r^[-](i;T) 



M-l 



-(T-t) 



{T~ty 



-$,,„[S](T-t). 



(47) 



The computational complexity to evaluate ^'|^^-'[S](i;T) for each i is 0(N^~^ ). Indeed, for each m, it is required to 
evaluate <I>i.m[S](T — t). This, in turn, requires at most N^^"^ evaluations of the Laplace transform of the symmetric 
Dirichlet distribution (one for each path (eo, . . . , e™) such that cq — e,; and Cj ^ Cj-i for j — 1, . . . ,m as seen in (4-1))- 



em) such that eg — e,; and Cj ^ ej. 
Therefore, the total computational complexity to evaluate (41) is 0{N^^^^) and, thus, the algorithm has complexity 



0{N^'') to compute the prices conditional on all starting regimes. We also remark that, for a fixed M , the computation 
o/4'i[S](0; T) can be sped up greatly by saving the paths (ei, . . . 
for each evaluation of $i^m[S](^) with C, & {5, . . . , k5}. 



) needed for (41) at the beginning and reusing them 



5.2 Error analysis 

We now proceed to give precise error bounds for our proposed method. First, we need to extend the framework of our 
algorithm. In this section we fix positive integers M and k, and an integer r e {1, . . . ,M}. We also let 6 = T/k and 
l(ei) = 1, for alH = 1, . . . , A^. Suppose the following conditions hold true: 

(1) The functionals 

$i,o[S](u;w),...,$j,r-i[S](u;w) 



defined in (26) can be computed exactly for any payoff function ^ : {ei, . . . ,eAf} — > M and any < u < v < T. 
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(2) For any m G {r, . . . , M — 1}, i G {1, . . . , N}, S : {ei, . . . , ejv} — > Hi, and < 6 < T, there is an approximation 
$i^„i[S](5; T) of <i>j^„j[S]((5; r), linear in S, such that 



$,,„[S](J;T) - $,,„[S](J;r) < Bmax|S(e,)| for ah 6 < So, 
$,,™[S](5;T)| < $,,™[l](5;T)max|S(e,)|, for all ,5 < Jq, 
for some < Sq < T and j > 1 and a constant B < oo, independent of S, i, m, and S. 
Next, define the following approximation for the pre-default payoff function 'i>i[E]{t;T): 



(48) 
(49) 



m=0 



ml 



-$,,,„[S](r-i;T). (50) 



Note that the case r = M essentially amounts to approximating the infinite series approximation (25 1 by the summation 
of the first M terms and, hence, that the M terms <i>i_o[S](^^; ^^)j ■ • ■ i ^i,M~i[^]{'u; v) are computed exactly. 

Consider the following generalized algorithm to approximate ^'[S](0;r) := (^E'l [S](0; T), . . . , \E'Ar[S](0; T))', which 
extends our approach in Section [O] 



(A) Approximate ^'j[S]((fc— 1)(5; T) by the functional vj/^^^'') [S]((fc— 1)(5; T) in (50 ); denote this first step approximation 

$(M,r.,fc)[2]((fc„ 1)J;T); 

(B) For each j from fc — 2 to (going backwards), compute iteratively 

^ $(M.'-)[sO-)](j5; {j + 1)6), (51) 

where S*--'^ is computed as 

The following result will be important in the subsequent error analysis. Its proof is given in the Appendix [A} 

Theorem 5.5. Under the above conditions (l)-(2), the previous algorithm (A)-(B) will result in a price approximation 
$(A^^'-,'=)[2](0;r) such that 



max 
ie{i,...,N} 



vI/aS](0;T)-vl,(*^''-''=)[S](0;T) 



< D- 



(52) 



where a := (r + j) A M and D — ||S||oo(l + B) with B defined as in (4-8) and ||S||oo '■— max£=i....^Af |S(£ 



The following result relaxes the condition (49), but does not yield a bound as explicit as in (52 1. 



Corollary 5.6. Suppose that the above conditions (l)-(2) hold true, but replacing (4-9) with the following milder 
condition: 



$,,™[S](^;T) < Cmax|S(ef)|, for all 5 < 5^, 



(53) 



for some < Sq < T and a constant C < oo, independent of 6, i, m, and :^. Then, there exists a constant E such that, 
for all k > 1, 



max 

i6{l,...,Af} 



< E- 



vl/,[5](0;T)-v|/f^''-^'=»[S](0;T) 

lfc-ti,...,iV> 

where a := (r + j) A M . 

We now consider different applications of our previous general approximation framework. 



(54) 



12 



5.2.1 Error analysis when the $i.m's are exact 



As a direct consequence of Tiieoreni |5.5[ by taking r = M and 5 = 0, the approximation 

M-l 

' " -$,;,™[S](T-t), 



m=0 



attains the following error: 



As indicated in the Remark 



max 

i6{l,...,Ar} 



< 5 



5.2 



for time-invariant generators, ^i^oi-^] ^^nd <i>i i[S] can be computed in closed form. 
Hence, using only these two values, the bond prices can be evaluated up to an error of order T^/k using k iterations 
and with a maximal polynomial complexity of 0{N^). Furthermore, even in the case of time-dependent generators. 



the functionals ^ifl[^]{C',T) and T) can be computed by the expressions (27) (at the cost of computing a 



single integral) and, thus, the same error will hold true. This last observation partly justifies our claim in Remark 5.3 



5.2.2 Error analysis for smooth time-dependent generators 

In the case of time-dependent generators, a natural approach is to approximate $i.m[2]((5; T) when S is small by freezing 
the values of f and K at the beginning of the integration interval [T — 6,T]. Concretely, consider the approximation 

"1 rSU,, 



$,,„,[S](<5;T) :=E^ 



-"(1+1) 



S(X„0 exp <j - E r{T - 5)' X^ds - ^ Xl^_^K{T - 5)Xn 



(55) 



I, n=0''"'^(") n=l 

These approximations can be expressed in closed form in terms of "symmetric" Dirichlet Laplace transforms as in 



Proposition 5.1 Obviously, (55) satisfies (53). In order to verify the validity of (48), let us assume that for all p, g, the 



functions fp and Kp^q are continuously differentiable on [0,T]. Denoting the random variable inside the expectation 



on the right-hand side of ( |55[ ) by G and from the definitions \2Q\ and (55), note that 
$,^„jS](<5;r)-$,,,„[S](<5;T) 



<S||S|U<5, 



for some constant B independent of S. Thus, (48) is satisfied with j = 1. In particular, if we use the approximation 



(55 1 for 771 = 1. . . . , Af — 1 (so that r 1), then the parameter a of Theorem 5.5 will be two and the resulting pre-default 



claim price approximation will attain an error bound of the form: 



max 
ze{i,...,Af} 



< E- 



5.2.3 Error analysis arising from the Taylor approximations 

Another natural issue is to quantify and control the error arising from using a Taylor approximation to the Dirichlet 



EH 



Laplace transforms (40) as suggested in Remark 5.2 Concretely, let us recall the following expression from Proposition 

1 



<i>.,™[s](C) 



[N - 1)™ 

(<!l,--,>!m) 



S(e„)e" 



■Qr Gt, 



eo 



! (^m— 1 Cm)); 



(56) 
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valid for time-invariant generators. Next, we consider the following "Taylor approximation": 

1 



(N-l) 



— J2 mra)e'^'''-~^--^ '^"-^'^"£,„(Cf'(go - -era), ■ • • , C^~'(g™-1 " Sm)), (57) 



where we approximate the Laplace transform (40 1 by the following Taylor polynomial of p -order: 



^ f 



AY 



dX. 



Let us now verify the conditions for Corollary 5.6 It is easy to verify that <i>i_m[^] satisfies the condition (531 . It 

ilities, 

X\\'dX 



remains to show the validity of (48 1. First note that, by the Cauchy-Schwarz's and Jensen's inequalities, 

i(x,A)r 



|/:,„(xi, . . . - . . . < m! / 



i=p+l ' 



< ml m 



i/2-l 



i=p+l 
oo 

— 5Z ™ 

i=p+l 



i=p+l 



where we had also used the fact that m! Jj, X\dX coincides with i*'' moment of a Beta(l,TO) random variable, which 
is known to be mli\/{m + j)!. Using the previous estimate, it follows that 



for an absolute constant K. Finally, its is clear that 



(58) 



(ei.....e™) 



Crr 



1 



(N-l) 



(=1..-..=™) 



can be bounded by ||S||tx)i?, for some constant R. In light of (58 1, Corollary 5.6 implies that the previous Taylor-based 
approximations will attain an error of the form 



max 
ie{i,...,iv} 



^,[S](0;T)-*|^'"^'=^^[S](0;T) 



< E- 



(59) 



with a := r + (p + 1) = (3 + p) A M if we take M > 2 and use the exact values (42 1 for m e {0, 1} and Taylor 
approximations of p*'' order for m = 2, . . . , M (so that r = 2). 



5.3 Pricing of decomposable claims on the underlying Markov process 

We proceed to extend our approach above to price a class of path-dependent claims, termed self-decomposable claims, 
whose payoffs can be decomposed into shorter maturity payoffs. We illustrate our method for a type of barrier option 
on the Markov process {Xt}t>o- As it is evident from (45), our approach heavily relies on being able to decompose 
the payoff of the claim into payoffs of shorter maturity. The following broad definition attempts to give a more precise 
meaning to this concept: 
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Definition 5.1. Consider a family of payoffs {^t.T}o<t<T, where for each < t < T, T^i.t represents a payoff 
depending on the path of X on [t,T]. We write T,t.T ■= ^i{Xs}t<s<T)- We say that the family {St.T}o<t<T *s 
self-decomposable if, for any < t < t' < T, the following decomposition holds true: 



^t,T — /(St 



for some measurable functions f,g'. 



We now proceed to describe our method. Following our strategy for simple claims, it is natural that a feasible 
procedure to price a self-decomposable claim Tit^riX) = T,{{Xs}t<s<T) will consist of the following two general steps: 

(Decomposition) Fix S = {T — t)/k for a positive integer k and apply the following decomposition with t' := t + 6: 



X, 



X, 



EQ [e^M'/(Et,tOS(XtO|^t =e,] + E^ e^M'g(E,_,,)S(XtO 



Xt 



where s(e;) := E^[e^*'-^ |Xf/ = e;] and s(e;) := E^^ [e^<:' •^'Ef ^t\ Xf = ei] . We then repeat the above decomposi- 
tion to evaluate the payoffs S(-) until T — t' \s small enough. 

(Near-expiration approximation) We proceed to apply an efficient approximation to evaluate claims of the form 
EQ [e^t-f h{Y.t^t')='{Xt')\ Xt = e,], when t is close to t' . 

Remark 5.7. The above method can also be extended to deal with claims whose payoffs can be decomposed in terms 
of the payoffs of other types of claims. For instance, we can say that two families of payoffs, say {s|"^}o<t<r and 
{Ej^^}o<t<T; o-fe mutually self-decomposable if, for any t < t' < T, 

for each k = 0,1, and some measurable functions fk, gk,hk,ik- 

As an illustration, we now consider the risk-neutral pricing of European barrier and digital contracts written on 
the volatility process {(Tt)t- One may view the process at '■= u' Xg as a proxy to the volatility of a market index. 
Instruments written on this process may be used to hedge volatility risk associated with periods of macro-economic 
bust or boom akin to that experienced by the U.S. economy leading into the 2008 crisis. Let us define the following 
family of path-dependent payoffs: 



I max a'Xs > Sj' 

^t<s<T ' 



and 



I max a'Xs < B]' 

''t<s<T ' 



The following simple relationships show that the family of payoffs {St,T}o<i<T and {St,T}o<i<T are self-decomposable: 

{1) I]t,T = ^t,t'^t',T, {2) j:t,T = Et,t'+Et',T-Et.t'Et',T = Et,t' + {l-Et.t')Et',T, {t<t'<T). (60) 

Now, let us consider the following European knock-out style barrier option: 



*f^^[S](t;r) :=EQ fe-^'^'^^''^St,TS(XT) 



Xt 



E^ 



Xt 



where we had used the same change of probability measure Q as in (23) and the following process analog to (44): 

/„_„:=-/ r'Xsds- ^ X'^-KX,. 
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Then, the foUowing decomposition follows from (60|: 



Xt = ei 



Xt = e. 



(61) 



For a near-expiration approximation method for '^f^^\p]{t;t') (i.e. when ^ := t' — t w 0), we use again (25), (26), and 



(411 



M-l 



OUT!'-: 



m— 



with 



1 



(TV - 1)' 



em)l{maXj(£T'ej)<B}e ^" = i ''"'^ £,„ (C?"' (§0 " e„), . . . X^' - em)). 



exp |— Jj"^ ^sdsl Et,TS(XT)|Xt = Ci can be easily 



Note that the knock-in style barrier option *f^[S](t;T) E*^ 

computed using the relation Sj^t = 1 ~ 5]^ In general one only needs to price either a knock-in or a knock-out 
contract, as the value of the other follows immediately from the knock-in/knock-out parity. 

5.4 Pricing of vulnerable call/put options 

We now develop a new method to price vulnerable call/put options on the stock. The latter type of options have 
received growing interest in the literature (see the introduction for further references). We set the terminal payoff 



q{St) — {St — K)^ in Eq. (36), thus obtaining a vulnerable time t call option price: 

n,(i;s) := n,{t,T;s) = \ {St - X)+e- Xt = e,,St = s 



As before, we will change the probability measure into Q so that, in terms of the process If^r defined in (44 1, 

U,{t, T; s) = EQ [ {St - K)+e'*-^ \ = e,, St = s] . 



Under the risk-neutral measure defined by Eqs. (34) and ([35|, we have that the stock price is given by 



St = St exp <^ J b'X.ds 



a'XAW? 



Elliott et al. 



(20051). 



where the evolution of Xt is determined by the risk-neutral generator Ay- under Q (see, e.g. 
Above, b := (6i, . . . ,1)^)' is given by hi := — erf /2. In particular, given (t{Xu '■ t < u < T) and St = s, we can 
see {Su}t<u<T as a geometric Brownian motion with a deterministic time-varying volatility and initial value s. As 
it is well-known, one can express the call price for such a model in terms of the Black-Scholes formula with constant 
volatility ct, short-rate r, spot price sq, maturity and strike K: 



BS {C,so,a\r,K) 



Concretely, denoting 



1 - 1 ~ 

J, _ ^ j {a'Xu)'^du, b ^ _ ^ J b'X^du, sq := se 
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we have 

E^[{St - K)+\a{Xu ■.t<u<T),St^s]^E^ {se 
Then, we obtain 

U,{t,T;s) = [e^*.^ {St - K) + \(j{Xu ■■t<u<T)] 

= EQ [e'''^BSiC; so, a^O, K) \ Xt ^ e„ St ^ s] . 



-K)+ =BS(C;so,a2,0,K). 



Xt = Ci, St = s 



Let us now focus on the time-invariant case, where is time-invariant and, hence, IIi{t,T]s) — 11^(0, T — t;s). We 
set 

rC 

F(C;s,e,):=EQ[(5c-i^) + |5o-s,Xo = e,]=EQ[e^-<Bs(C;seJ■^''^"''^^' / {a' X^)^du,0, K] S'o = .s,Xo-e, 



Our approach is based on two ideas. Firstly, if C is smah, note that 



F(C; s, e,) « F(C; s, e,) := E^ [e^«'<BS (C; se^'''^^ , ((t'Xc)^, 0, i^) 



Sq = s,Xo = Ci 



(62) 



up to an error 0(C)- This is because /^^ r' Xudu — C,r' Xq and /p^(cr'X„)^(iM = C,{a' X^Y if there are no transitions of the 
process (Xt) during [0, C]. Since the expression in (62) can be seen as a European claim of the form (201 with maturity 
C and payoff S(ej) :— BS (C; se'-''^ , ct|, 0, i^T) , one can evaluate this "first order approximation" using our Algorithm [l] 
in Appendix [B] 

For a general maturity C, we proceed as in (45 1. Concretely, for (5 < C, we have the recursive relationship 



F(C; s, e,) = EQ [e''>-'F{C - S; Ss, Xs)\ So = s,Xo = e,] = E^ 



Xo = Ci 



where now (C — (5) ^ [a' X^)'^du and 6 := (C — <^) ^ J^. b'Xudu. As before, the last expression is approximated 

by 



i^(C;s,eO :=E^ 

= eQ 



^5 



d: se 



Sb'Xs + VSa'Xs 



-^0 = Si 

1 



/27r 



ATo = Ci 



In principle, for a fixed s, one can see the right-hand side above as the price of European claim of the form (20 1 with 
maturity S and payoff S(ej) := E^ F (^(^ — 6; se'^^^+°'j^^?, e^^ . But, since this approach would require to compute 
i^(C — S;p, Bj) for all p and e^, this would be computationally inefficient. To resolve this issue, we restrict all possible 
initial prices s to be in the lattice £a.b '■— {se^'^ : i G {—B, —B + 1, . . . , _B — 1, B}} for a small A and a positive integer 
B. Then, we can approximate F{(; se^^ ,ej) as follows: 



f(C;se''^,e,)«EQ[e^«'^S,(X5) 



Xo = Ci 



with 



where 



k=-B 



1 

k-1 ^ 



''^dz. 



y J . - 



{k - i + 1/2) A -Sbj 



{k^-B,-B + l,...,B~l), z%_^ 



Note that the points 2:^'-"s are chosen so that the midpoint z^'-* of the interval [z^'f.]^, z^:"*] is such that sexpjfcA} 
sexpjiA -|- 5hj + ^5i7jz]f} and, hence, above is a Riemann-Stieltjes sum approximation of the payoff 

1 



5,(ej) 



F 



27r 



'/^dz. 
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6 Numerical Simulations 



The objective of this section is to assess the accuracy and computational speed of the different methods introduced in 
the previous section. We first describe the simulation scenario in Section |6.1| and then provide the numerical results in 
Section EH 



6.1 Simulation Scenario 

Throughout this part, we use the following parameter specifications borrowed from the parameter setup of |Capponi| 



and Figueroa-Lopez (2012). The right panel also shows the loss rates corresponding to three distinct regimes, hereafter 
referred to as "low" , "middle" , and "high" default regime. The default intensities are based on the empirical result of 
( 201l[ p In agreement with empirical market evidence, we have selected the loss rates to be increasing 



Giesecke et al. 



in the credit riskiness of the regime. 





1 (low) 


2 (middle) 


3 (high) 


1 


-0.380313 


0.33687 


0.043443 


2 


0.254397 


-0.254397 





3 


0.208683 


0.000006 


-0.208689 





h 


L 


1 


0.741% 


10% 


2 


4.261% 


40% 


3 


11.137% 


90% 



Table 1: Left panel shows the generator of the Markov chain under the risk- neutral measure. The rows indicate the 
starting state, while the columns indicate the ending state of the chain. The right panel shows the default intensities 



as reported in Giesecke et al. (2011 1 as well as our loss rates given default associated to three regimes. 



6.2 Numerical Results 



We start by illustrating the accuracy and computational speed of the method described in Section |5.1| for computing 
the pre-default price of a defaultable bond. These bond prices, denoted hereafter ip{t) := (V'i(i)j • ■ • , i^N{t))' , are given 
by the formula (20) with payoff function S(ei) = 1. Concretely, we define 



(63) 



We compute ( 63 ) by our method and a standard numerical solution of the Feynman-Kac representation ( 29 ) for the 
bond price ( |63[ ) . Table [2] shows the time-0 bond prices for different maturities computed using our method and a 
Runge-Kutta type numerical solutiorj^ of the system ( 29 ) with terminal condition '^i (T) — 1 under the parameter 
setup of Table [I] Figure [I] shows the time t bond price for maturity T — 1 year (left panel) , and T = 20 years 
(right panel). It is evident that our method is highly accurate even for long maturity bonds. Furthermore, according 
to our computational experiments, our method is in most cases more efficient than either solving the ODE system 



by the Runge-Kutta algorithm or computing the exponential ( 33 ) using a Fade type approximation. For the sake of 



completeness, Appendix [C| compares the processor time for these three methods. 

Next, we illustrate the power of the algorithm to compute prices of the knock-out digital and call options described 
at the end of Section |5.3[ Table [3] reports the prices of knock-out digital and call options written on the volatility 
process obtained by applying our methodology, and compares them with the corresponding Monte-Carlo estimates. 
The results show that our method is highly accurate even for long maturity options. We conclude the numerical 
analysis with an illustration of the algorithm performance on pricing vulnerable European type claims, which depend 



• ^Giesecke et al.| | |201l| employed a three-state homogenous regime-switching model to examine the effects of an array of financial and 
macro-economic variables in explaining variations in the realized default rates of the U.S. corporate bond market over the course of 150 



years. 

^This solution was obtained using the MATLAB function ode45. 
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^0 


= ei 


^0 


= 62 


Xo 


= 63 


T (yrs.) 


NM 


ODE 


NM 


ODE 


NM 


ODE 


0.25 


0.9921 


0.9921 


0.9884 


0.9884 


0.9686 


0.9686 


0.50 


0.9837 


0.9837 


0.9772 


0.9772 


0.9393 


0.9393 


1.00 


0.9659 


0.9659 


0.9555 


0.9555 


0.8864 


0.8864 


2.00 


0.9282 


0.9281 


0.9146 


0.9146 


0.7991 


0.7990 


5.00 


0.8140 


0.8136 


0.8029 


0.8031 


0.6274 


0.6273 


10.0 


0.6488 


0.6484 


0.6430 


0.6431 


0.4701 


0.4701 


15.0 


0.5170 


0.5166 


0.5131 


0.5131 


0.3691 


0.3690 


20.0 


0.4119 


0.4116 


0.4090 


0.4090 


0.2931 


0.2930 


25.0 


0.3282 


0.3280 


0.3259 


0.3259 


0.2334 


0.2333 


30.0 


0.2616 


0.2613 


0.2597 


0.2597 


0.1859 


0.1858 


50.0 


0.1055 


0.1053 


0.1047 


0.1047 


0.0750 


0.0749 



Table 2: Time t — bond prices for different time-to-maturities T using the ODE method and the new method (NM) 
with parameters M — 2 and S — 2.5 years (see Algorithm [l] in Appendix [b|). 



on both the diffusion and the Markov chain component. Figure [2] compares our novel method with a plain Monte 
Carlo method, given that the initial regime is 1 (left panel), and 3 (right panel) (regime 2 is quite similar to regime 
3). As seen there, the new method significantly improves the quality of the approximation compared to the first 
order approximation (62), especially for longer maturities. For the sake of completeness, we have included the precise 
algorithm in Appendix [b] (see Algorithm [2] therein) . Error analysis and further extensions of this method will be 
postponed for a future publication. 



7 Conclusions 



We have rigorously analyzed pricing of vulnerable claims in regime-switching markets, and provided novel semimartin- 
gale representations . The main theoretical tool underlying our results is a novel Poisson series representation of a 
generic claim price, whose payoff depends on the underlying Markov chain. Using this representation, we were able 
to prove differentiability of the pre-default price function under mild conditions on the Markov chain generator. This 
allowed us to provide rigorous Feynman-Kac and semimartingale representations of the price process of the vulnerable 
claim. We then constructed a novel pricing methodology which makes use of the Poisson series representation, and 
expands the price of the vulnerable claim in terms of the Laplace transform of the symmetric Dirichlet distribution. 
We provided a detailed error analysis of our proposed algorithm under different sets of assumptions relating to how 
each conditional price component may be approximated. One of the novel features of our method and error analysis 
is that it also applies to time-dependent Markov-chain generators. We further demonstrated the robustness of our 
methodology by showing that it can be used to price path dependent claims, termed self-decomposable, as well as 
vulnerable European-type options, whose payoff depends on both the diffusion and the Markov-chain component. We 
numerically demonstrated the algorithm on defaultable bonds, barrier options on the stock volatility, as well as vul- 
nerable European call option prices. In all these cases, we have shown that it achieves a high level of accuracy, and, at 
the same time, it is computationally fast when compared to standard PDF or Monte-Carlo based methods. 
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Figure 1: Bond price comparison for a one- year and a 20-year bond using the ODE method and the new method 
(NM) parameters M = 2 and S = 2.5 years (see Algorithm [T] in Appendix [b|) . 




Figure 2: Comparison of vulnerable call option prices using Monte Carlo and the new method of the Algorithm [2] in 
Appendix [B] for maturities of T = 3 years and T — 6 years. Above, k = 2 refers to the number of discretization points 
on the PriceClaim Method of Appendix |b] (i.e. 5 ~T/k therein). 
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Knock-Out Di 


jital Contracts 


Knock-Out Barrier Call Options 






5(Xt) = 1 






5(Xt) := 


(fTx^ - K) 








B = 


= CT3 




B = a3, K ^ 0.075, Units = 


= 10 




^0 


= ei 




= 62 




= ei 




62 


T 


NM 


MC 


NM 


MC 


NM 


MC 


NM 


MC 


0.50 


0.9719 


0.9693 


0.9525 


0.9523 


3.5081 


3.5574 


21.1632 


21.0587 


1.0 


0.9419 


0.9455 


0.9093 


0.9100 


5.8869 


5.7910 


18.2870 


18.2697 


2.5 


0.8482 


0.8466 


0.7980 


0.7964 


9.0631 


9.0351 


13.1068 


13.0927 


5.0 


0.7013 


0.6986 


0.6506 


0.6481 


9.1533 


9.2635 


9.3528 


9.5317 


10.0 


0.4732 


0.4721 


0.4373 


0.4336 


6.4919 


6.3608 


6.0298 


5.9426 


15.0 


0.3187 


0.3201 


0.2944 


0.2955 


4.3833 


4.3335 


4.0508 


4.1597 


20.0 


0.2146 


0.2141 


0.1983 


0.1974 


2.9521 


2.9365 


2.7274 


2.8137 


25.0 


0.1445 


0.1445 


0.1335 


0.1341 


1.9879 


1.9577 


1.8366 


1.7792 


30.0 


0.0973 


0.0985 


0.0899 


0.0912 


1.3386 


1.3136 


1.2367 


1.2111 


35.0 


0.0655 


0.0676 


0.0605 


0.0608 


0.9014 


0.9038 


0.8328 


0.8228 



Table 3: Knock-out digital and call options on the volatility process using the new method (NM) and the Monte Carlo 
method (MC). The call option prices are expressed on the 10"^ scale and the results for Xq = have been omitted 
since they are knocked out at contract initiation. Here, r = (0.01,0.1,0.3)' and a — (0.05,0.1,0.2)'. 
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A Proofs 



Proof of Lemma 12.21 

By the definition of Q and the fact that 77 is a G-martingale (see the paragraph before (15)), it suffices to prove that 
rjt^'^ is a G-martingale under P. From Ito's formula and the definition of in ( 15 1, we have the process 



Jo Jo 



From (|2|, (15), and (19), (»7tCf )t can be written as 



•'^ 0<s<tk.l=l 



where H^'^ :— X]o<s<t ^{x _=fe}l{Xs=i}- Since the first two terms on the right-hand side of the previous equality are 
(local) martingales under P, it remains to show that the last term vanishes. But, given that AHg 7^ at s = r, the 
summation in the last term above will be provided that AXr = 0, a.s. In order to show this, let us recall that by 
definition X has no fixed-jump times; i.e. V{AXt ^ 0) = for any fixed time t > 0. Also, using the definition of r 
given m Eq. ([l|, T = mi{t > : /J h{Xs)ds > 9}, where 9 is an exponential random variable independent of X. Then, 
conditioning on X, ¥{AXt 7^ 0) = E |E [ 1{aav#o} | s ^ O]}. Denoting < ri < T2 < . . . the transition times of 
the Markov chain X, E [1{aXt^o}| s > O] is given by 



E 



=n} 



Xs,s>0 



= E 



El 



X,,s>0 



= 0, 
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where the last equality follows from the independence of X and Q, and the fact that is a continuous random 
variable. □ 

Proof of Theorem [3TT1 



For simplicity we omit the dependence on the claim's payoff ^ and maturity T in (20) and just write Recall 
that ^{t) (^'i(t), . . . , *Af(i)) is such that 



7(t;T) ^E'^ E{XT)e-'^t'^'''-+''^^''^'^'\Xt 



(64) 



But also, changing into the probability measure Q, we can write 



7(t;T) =£^2 \E{XT)e--^''^'''+''^^'^'^' 
and, hence, we have the following representation for 



Vt 



Xt = e,; 



(65) 



The solution of (21) can be written as 



Vt ■= e 



see pp. 334-335 in 



Bielecki et al. 



( 2008-b[ ). Let K{t) = [K,j{t)]ij and P{t) := (fi(i), . . . ,fAr(i))' be defined by 

N 

K,^j{t) := log (1 + K,j(i)) l^^j = - log ((iV - l)afj{t)'^ l.^j, f,(i) := r, + h,L,~ ^ a^^.(t)S,j(t). 



Then, since by the definition (14 1, AHl'^ = 1 whenever there is a jump from the state to the state Cj (i.e. Xg 



and Xg = Cj) and AHl'^ — 0, otherwise, we have 

^log('l + ^^,^,(s)Ai7:'A =^^log(l + ^.,,(s))Aif:-^= K-K{s)Xg 



Similarly, by the definition (14), 



2J 

where v = {vi, . . . , vn)' with Vi{s) := X^^i Then, since Vg + hgLg = f'Xg with f = (ri + . . . , tat + 



HnLn), we obtain from (651 that 



exp ■ 



where X' denotes the transpose of X. 



?{syxgds- K-K{s)Xg 



Xt = Ci 



□ 



Proof of Lemma 14.11 

Again, for simplicity we omit the dependence on the claim's payoff 2 and maturity T in the functions \E'i[S](i;T) 



and <i>i.m[S](t; r). In light of the series representation (25), we see that it suffices to show the function $i.m(C) is 



continuously differentiable in C G (0,r) for each m > 0, and that there exists a sequence {Km}m>o of reals such that 



(i) sup |$,„(C)| < (ii) sup K,(C)|<ifm, (iii) V e-^^/C™ <cx). 

0<C<T 0<C<T — 



m=0 



22 



We now show that (i)-(iii) are satisfied provided the function K and f defined in (24) are such that 

sup |ri(s)| < oo, sup \Ki_j{s)\ < CO, and sup jif,- j(s)| < oo. 
s6[0,T] se[o,T] ' se[0,T] 



(66) 



The conditions (66 1 directly fohow from (22) and (28 1. Denoting Vm,i{C) the random function inside the expectation 

(67) 



in ( 26 ) , one can check that 



sup |7'™,^(C)| < S(ei) • e"^'"^'^'""Po<C<T 1*^.(0'!+™ max.,, supo<j<r|if.,j(C)l = 2(60^"", 
0<C<T 



for a constant A < oo. Also, Vm.iiC) is continuously differentiable and 

{m m 
J2 [U(n+i)r{T - Cf>(n+i))'^n - t/(„)r(T - CC/(„))'^n] - ^ K-iK' {T - CC/(„))X„ 
n=0 n=l 

where := 1 — ?/(„). In particular, there exists a constant B < oo such that 

sup |n'„(C)l<Sm^". 

0<C<T 



in) 



(68) 



By the formal definition of the derivative $^i(C) ^ind the Dominated Convergence Theorem, one can check that 

□ 



(|67|-(l68]) will suffice for (i)-(in). 



Proof of Theorem l4T2l 

The following three parts show each of the statements of Theorem |4.2[ 
(1) Consider the (F,Q)-local martingale 

Z{t) :=EQ \E{XT)e^-fo^^''^+''^^'^'^'\Tt 

Note that 



and, thus, by Ito's formula and the decomposition (16), we get 



dt 



(69) 



where we have used the fact that ^'(t) is differentiable in t by virtue of Lemma 4.1 Since Z{t) is a (F, Q)-(local) 
martingale, its drift term must be zero. When Xt — Ci, such condition translates into the coupled system of differential 
equations for the price function given in Eq. (29). 



(2) The following dynamics for Z also follows from (69) and the martingale property of {Zt} 

dZ{t) ^ e--/'o('-"+''"^")''"(*(t),dM<2(i)). 
In particular, the price process 7(t) follows the dynamics 

d-i{t)^d(e^>''^^^^^^'^''Z{t)^ 

^ {n + htLt)e^o(^^+h^i^'')dsz(^t)dt + efo(^s+KL^)ds^2;{t) 
^ in + htLt)e^>-'+''^^'^'^''Z{t)dt + {m{i),dM^{t)) . 



t>o- 



(70) 
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Next, denoting H{t) := 1 — H{t) = lr>t and applying Ito's formula, 

dr{t; T) = d{H{t)j{t)) = H{t-)dj{t) + j{t-)dH{t) + A'y{t)AH{t). (71) 
By the continuity of ^'(t), 

A^it)AHt = (*(t), AXt) AHt = (*(t), AXr) = 0, a.s., 

where the last equality hold true in light of the construction ([!]) and the independence of X and the exponential variable 
8 there (see also the proof of Lemma 2.2 above). Combining ([70| and (71 ), 



dT{t;T) = i/(f-)(rt + /itLOe-''«^''=+''^-^=^'^'^(i)dt + -ff(t")(*(0,f^^"^W) +7(i")rf^W 
= {rt + htLt)T{t; T)dt + H{r) (t) , dM'^ (t)) - -f{r)dH{t) 



By the Doob-Meyer decomposition of Lemma 2.2 we have dH(t) — dQ^ + (1 ^ H{t ))htdt and, thus, 

dV{t-T) = (h + htLt]T{t- T) - -i{t-)ht[l - H{t-)]) dt + H{t-) {'f{t),dM'i{t)) - j{r)d^f 

The last step follows from the fact that on the event t > t, we have T{t'; T) — ^{t^\T), by definition of r(t, T), and 
also 7(t~) — {'^{t), Xf-) . Therefore, we can write the risk-neutral dynamics of prices process T{t;T) as 



- {n + ht{Lt - l))dt + ^^^^^ ^^^^ -dg, t<T, 



nt';T) 

which concludes the second statement of the theorem 



(72) 



(3) The price dynamics of the vulnerable claim under the real- world measure follows directly from plugging Eq. (17) 
into ( 72 1 to get: 

dr{t,T) 



T{t-,T) 



= rt + ht{Lt- 



^ {^it),iA'it)-iA^y{t))X,) \ ^ {<f{t),dM^it)) 



*(i),^t-) 



where we also used that ^ = (see Lemma 2.2). Finally, we get the dynamics (31 1 since 

<vE'(t),(A'(t) - [A^)\t))X,) {{A(t)-A^{t))m.^i) 



When Xt — et, the previous quantity simplifies as 



A(t) :=E(«^^ W-«?.W) 



*(^),^t> 
Ml 



(73) 



and ( 73 ) can be written as ((Z?i(t), . . . , 0^(1))' , Xt). 



□ 



Proof of Theorem 14.4 



The following three parts show each of the statements of Theorem 4.4 
(1) Consider the (F,Q)-local martingale 



Note that 
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Thus, by Ito's formula, we get 



where we used the notation -f{t; St) = g{ST)e- -f* (r.+/i.L,)ds|jr^ 

dZ{t, St) = -{n + htLt)e- fo(r^+h.L^)dn (n(t; St) , Xt) dt + (n(t; St),Xt) . 

Under the stated differentiabiUty assumptions, we have that 

d{n{t-St),Xt) ^ C{U{t;St),Xt)dt + dMt, (74) 

where C is the generator of {t,St,C't), and {Mt}t>o is the local martingale component of {t^St,Ct), under the risk 
neutral measure Q. These have been given in a more general context in Capponi and Figueroa-Lopez (2012), Appendix 
A. In our specific setting, they reduce to 



LVVi(t; s) = — h r.iS- 



dt 



ds 



52 s 



and 



N 



dMt^Y.{Y. ^*-) - St-)] dMl'' + a^St ^^'^^^^''^ HldW^ 

i=l ^ j^i 



(75) 



(76) 



where HI and have been defined in Eqs. (14) and (13), respectively. Using the operator £ and local martingale 
Mt given by Eq. (75) and (76), respectively, we obtain 



dZ{t, St) = -{rt + htLt)e- Io(^^+'^^L^)du ^^(t; St),Xt) dt + Io(r..+h.L^)du ^^(t; St), Xt))dt + dMt] 



( j_ w wmf Q ^ y \ ^ f^ffll^^ll^ ^ ^dm;St),Xt) 
(n + htLt) {U{t; St), Xt) H 7^, ^ nSt- 



dt 



ds 



ajS^d^ m-St),Xt) 



- a%^^{t){n,{t-St)~nc,At-St)) 



dt + e-ioi^-^h^^-'^du^Mt. (77) 



2 d^s 

As {Z{t, St)}t>Q is a (F,Q)-local martingale, its drift must vanish. When Xt — Ci, this translates into the coupled 
system of partial differential equations given by Eq. ( 36 1 . 

(2) In light of the martingale property of Z{t, St), we have that its dynamics must be 

dZ{t,St) = e-J'o(''"+'*-^")*'dXt. 



Using results in Chapter 11 of Bielecki and Rutkowski (2001 ), we can re-write the first term of dA4t as 



AT 



^ J2 [n,(i; St) - n,(i; St)] dMl^' = (n(i; St),dM^), 
1=1 j^i 



(78) 



Moreover, from the definition of HI — l{Xt=ei}j given in Eq. (14), we have that the second term of djvit can be written 



as 



N 



dMt; St) 
ds 



HidW^ = (TtSt 



d{U{t;St),Xt) 

ds * ' 



and, thus, altogether 



dZit, St) = e-/oK+''"^")*' ( /n(t; St), dMf \ + atS, 



d{U{t;St),Xt) 
■ ds 



dW> 



(79) 



(80) 



In particular, the process 'y{t;St) — e^ai^^+^^^^^)'^^ Z{t, St) follows the dynamics 
d^{t-St) = (r, + /i,L0e/o(-+''"^")'^"Z(<,50di+ (n(i;5,),dM,«\ 



ds 



dW} 
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Denoting H{t) := 1 — H{t) = lr>t, we have that the vulnerable claim price process r{t,T;St) := H{t)j{t; St) has 
dynamics 

dTit, T- St) - d{H{t)^{t; St)) = Hit-)dj{t; St) + i{t- St)dH{t) + l\^{t- St)^H{t). (81) 

where, using the hypothesis that n(t; St) is continuously differentiable in t, and twice continuously difFerentiable in St, 
we have 

A7(t; St)^Ht = (n(i; St),^Xt) AHt - (n(t; St), AX,) = 0, a.s. 

with the last equality being true in light of the construction ([ij and the independence of X and the exponential variable 
8. By the Doob-Meyer decomposition of Lemma 2.2 we have dH{t) = d^^ + (1 — H{t^))htdt and, thus, 



dTit, T; St) = ((rt + htLt)rit, T; St) - i{t; St)htH{t-)) dt + H{t-)atSt 
+ H{t) (n{t; St),dM^) ~ -i{t; St)d^f 



dm;St),X, 



ds 



'-^dW^ 



(82) 



On the event t > t, we have r{t,T;St) = l{t;St) by definition of r{t,T;St), and also j{t;St) = {n{t; St), Xt-) . 
Therefore, we can write the risk-neutral dynamics of prices process T{t,T; St) as 



T{t,T;St) ' '^'^ " ^" ^ ' {^{t;St),Xt 

which concludes the second statement of the theorem. 



m\St),Xt-) 



(3) The price dynamics of the vulnerable claim under the real- world measure follows directly from plugging Eq. (17) 
into (831, and using the P Brownian measure defined in Eq. (35) to get: 

dV(i,T-St) 



T{t-,T:St-) 



rt + htiLt~l) + 



{n{t;St),iA'{t)-{A^y{t))Xt) , St 



d{n(t;St),Xt 



ds 



(Il{t;St),Xt 



{U{t;St),Xt^) 

(U{t;St),dM^{t)) 9imf),x.) 

' '^tSt ,^,^_ ^\ ^ , dWt-dil t<T, 



{fit - rt) ]dt 



{nit;St),Xt 



{Uit-St),Xt-) 



where we also used that = ^'^ (see Lemma 2.2). Finally, we get the dynamics (38) since 

{U{t;St),{A'{t)-{A^y{t))Xt) _ {{A{t)~A'i{t))U{t;St),Xt) 



{n{t;St),Xt) 

When Xt = Ci, the previous quantity simplifies as 



{nit;St),Xt 



(84) 



(85) 



D.{t,St) :=5](a,,(0-a^,(t))^5j||] 



and (85) can be written as {{Di{t,St), . . . , D^^t, St))' , Xt) . 



□ 



Proof of Proposition |5.1[ 

From the assumed time-invariance of the parameters, we have the following: 

$,,„[S](C) =Ef [s(l„,)e-5:"^oC(c/(„..,-r/<„,).'x„-E"..<-.i?x„ 

Let A„ :— — U(^n-i) foi' ^ = 1, . . . , m -I- 1 and note that X^ni^ ^ ^- well-known 
that the distribution of 

A := (Ai,...,A™) 



Kendall and Moran 



(1963) 
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is the symmetric Dirichlet distribution with parameter a := (1,...,1) G M"'+^. We recall that a random vector 
A := (Ai, . . . , A„i) follows a Dirichlet distribution with parameters a :— (ai, . . . , am+i) G M™+^ such that min^ > 
if its density is given by 

m 

D{X,, ...,Xm)= B{a)-' Yl - Ai A™)"'"+^-11t„(Ai, . . . , A„), 

1=1 

where B(a) := UZV {E"=i^ and T„ {(Ai, . . . , A,„) G : A, > & Y:T=i ^^ < !}■ Next, condition- 

ing on {Xi, . . . , Xjn) and using independence between U and X, 

$,,„[S](C) = Ef [ / e- ^"^0 ^^'^"^"+^i?(Ai, . . . , A™)dAi . . . dX^E{X„,)e- ^-^i 



where Am+i = 1 ^ X]™=i '^n- The following simplification can be made: 



e" E"=o C^'-*"^"+ii:)(Ai, . . . , A„)dAi . . . dA„ 



where Cm is given as in (40 1. Finally, note that, by construction, 
and, thus, 



E-o^ C.'(x„-x„OA„..-f.'x,„^(^^^ ^^^^^^ 

= efc|Xi_i = e„i) = 1/{N ~ 1) for any k ^ m 

ei, • • • = ejn\Xo = Ej) = _ j^^m ' 

for all (ei, . . . , e^) G {ei, . . . , e^v}™ such that ^ e,;_i. In that case, it is clear that 

(=1,...,^™) 

□ 



1 



Proof of Theorem 15.51 

Let l(ei) = 1 for all i = 1, . . . , A^, and ||S||oo — niax^ |S(e^)|. Let us first note that for all m < t; with v — u < S, 



-<i>,^m[l]{v-u;v) < \\E\\^^,[l]{v-u;v) < (86) 



m=0 



where for the second inequality above we have use (49) and the fact that < 1 (being the bond price under the 

z*'' economic regime). In particular, from the definition (51 1, 



max 



< max 



max 



< ■ ■ ■ < max 



^''^'■■''^mk^mks) <\\E\\ 



Thus, by backward induction, it follows that 



max max 

j=Q,...,K-l i=l,...,N 



< lisi 



(87) 
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Next, note that, for j — k — 1 and D ||S||oo(l + B), 



(M.k.r) r_,- 



S] U6; T) - [S] {36;T) = | *f ^''■^ [S] (j<5; T) - M', [S] (j<5; T) 

M-l „ 



Af-1 



m—r 



■ E 

7n—M 



where we used the representation ( 25 ) , ( 48 1 , and the fact that 

m=M ■ £=0 



Next, for j = k — 2, using (46) and (51 1, 



^(M,fc,r)j^j^^.^. _ vI/JsKjJ; T) < vl/(*^^'-)[s(^)](j5; (j + 1)5) - vl/4s(^")](j<5; 0' + 1)6) 



Ti 



(89) 
(90) 



T2 



where 



S(^)(e,) :=*f^^'^''^)[S]((j + l)<5;T), S(^)(e,) := vI;4S]((j- + 1)<5; T). 



The term in (89 ) can be bounded analogously as in (88 1, but replacing the constant D by ||^(-'' ||oo(l + S), which itself 
can be bounded again by D = ||S||oo(l + B) in light of (87). Thus, 

Similarly, we have 



To = 



^^[SO-) - S(j")](j<5; (j + 1)6) < max S(^')(e,) - S(^)(e,) 4'.[l](j(5; (j + 1)<5) < max E'-J\ee) - E<^^\ei) 



which can be bounded by in light of (S. Putting together the previous two estimates, 

^M,k,r) _ 2)<5; T) - *,[S]((fc - 2)6; T) 

Proceeding by induction, it follows that 



^p'=''")[S](M;r)-«',[S](M;T) < {k ~ 1)06^'''+^^''^'' = {k - e)D 



T 



{r+j)AM 



and, in particular, we obtain (52) by setting ^ = 0. 
Proof of Corollary |5.6[ 



□ 
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The proof is very similar to that of Theorem 5.5 The key difference is in the estimate (86 1, which needs to be modified. 
Under the modified condition (53 1, we proceed as follows, fixing = v — u, 

r-l 



(M,r) r^- 



m=0 



m— m—r 

< \\E\\^^^,m;v) + \\E\\^CC 
<\\E\\^{1 + CC). 



M-l 



In that case, by backward induction, Eq.(87l will take the form: 



max max 

j=Q,...,K-l i=l,...,N 



cr\fc 



(91) 



for a constant 1 < < oo (depending only on C, T, and r). Following the same proof as in Theorem 5.5 we will have, 
by induction, that 



«'f^'*^'''^[S](£<5;T) - *,[S](M; T)\ < {k ~ l)E5^^+'^''^' = {k - l)E 



{r+j)AM 



with E := ||S||oo(S + F) and Eq. ^ follows. 



□ 



B Pseudo-codes 



Algorithm 1 [Prices] = PriceClaim(r, ^, 6, M) 
C = mm{T, 6} 
for i = 1, . . . ,iV 

Prices(z, 1) = e-C E™=o S *.,™[S](C) 
endfor 
ii T <S 

Return Prices(:,l) 
else 

for j = 2,3,...,rfl 
for i = 1, ... ,7V 

Prices(z,.?) - e-^E™to £ 'i>Mn[Prices(:, j - 1)]((5) 
endfor 
endfor 
endif 

Return Prices(:, :) 



C Computational speed comparisons 

Table H compares the computational times for our proposed method (Algorithm [T]) with respect to both the Runge- 
Kutta and a matrix exponential approximation method based on the representation ( 33 1 using the parameter setup of 
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Algorithm 2 [Price] = ConiputcVulnOption(r, S, A,M,k) 
ii T <S 

for i = -B,-B + 1,...,B 
for j = l,2,...,iV,...,B 

endfor 

[Price(i, 1), . . . , Price(i, N)] = PriceClaim(T, E^,T/k, M) 
endfor 
else 

Price„eM) = ComputeVulnOption(T — 5, 5, A, M, k) 
for i = -B, -B + l,...,B 
for j = l,2,...,iV,...,B 

endfor 

[Price(i, 1), . . . , Price(j, N)] = PriceClaim((5, %,5/k, M) 
endfor 
endif 



Table [T] given in Section [63] For each given pair ((5, T), the bond prices '^p2{C,)^ "03(0) ^'I'g computed for ah time- 

to-maturity C, ^ {6,26, ... ,k5 := T}. The entries in Table |4] represent the ratio of the new algorithm's computational 
time with M = 2 to the respective ODE Method computational time and the matrix exponential computational time. 
As we can see, our method tends to outperform ODE Method for maturities shorter than 15 years and also for mesh 
sizes larger than 0.05 years; however, even in the worst case our method operates at a speed which is very comparable 
to the faster of the two. 
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Table 4: The top panel shows the ratio of the processing time using our proposed method (Algorithm [T] above) and 
a Runge-Kutta Type numerical solution of ( 29 1 using the MATLAB function "ode45" . The bottom panel shows the 



processing time ratio between our method and a Pade type approximation of the matrix exponential ( 33 1 using the 
MATLAB command "expm". 
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